Origin of asymmetric reversal modes in ferromagnetic/antiferromagnetic multilayers 
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Experimentally an asymmetry of the reversal modes has been found in certain exchange bias 
systems. From a numerical investigation of the domain state model evidence is gained that this 
effect depends on the angle between the easy axis of the antiferromagnet and the applied magnetic 
field. Depending on this angle the ferromagnet reverses either symmetrically, e. g. by a coherent 
rotation on both sides of the loop, or the reversal is asymmetric with a non uniform reversal mode 
for the ascending branch, which may even yield a zero perpendicular magnetization. 
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PACS numbers: 75.70.Cn, 75.40.Mg, 75.50.Lk, 85.70-w 

For compound materials consisting of a ferromagnet 
(FM) in contact with an antiferromagnet (AFM) a shift 
of the hysteresis loop along the magnetic field axis can oc- 
cur which is called exchange bias (EB). Often, this shift 
is observed after cooling the entire system in an exter- 
nal magnetic field below the Neel temperature Tn of the 
AFM (for a review see Q)- The role of the AFM is to 
provide at the interface a net magnetization which is sta- 
ble during reversal of the FM consequently shifting the 
hysteresis loop. The key for understanding EB is to un- 
derstand this stability. 

In the approach of Malozemoff ^ due to interface 
roughness domain walls in the AFM perpendicular to the 
FM/AFM interface are supposed to occur during cooling 
in the presence of the magnetized FM resulting in a small 
net magnetization at the FM/AFM interface. However, 
the formation of such walls and the stability of the inter- 
face magnetization has never been proven. Additionally, 
the formation of domain walls exclusively due to inter- 
face roughness is energetically unfavorable and therefore 
unlikely to occur. Therefore, other approaches have been 
developed where a domain wall forms in the AFM par- 
allel to the interface while the magnetization of the FM 
rotates 0,0. However, it was shown by Schulthess and 
Butler [5j that in this model EB vanishes if the motion 
of the spins in the AFM is not restricted to a plane par- 
allel to the film as was done in Koon's work. To obtain 
EB Schulthess and Butler assumed uncompensated AFM 
spins at the interface. But their occurrence and stability 
during a hysteresis loop is not explained, neither in their 
model nor in other similar models 0, Q . 

In a recent experiment Miltenyi et al. Q showed that 
it is possible to strongly influence EB in Co/CoO bi- 
layers by diluting the AFM CoO layer, i. e. by insert- 
ing non-magnetic substitutions (Coi_a;Mga;0) or defects 
(Coi_yO) not at the FM/AFM interface, but rather 
throughout the volume part of the AFM. In the same let- 
ter it was shown that a corresponding theoretical model, 
the domain state (DS) model, investigated by Monte 
Carlo simulations shows a behavior very similar to the 
experimental results. According to these findings the ob- 
served EB has its origin in a DS in the AFM which occurs 
during cooling and which is stabilized by non-magnetic 
defects carrying a net magnetization at the AFM inter- 



face. Later it was shown that a variety of experimental 
facts associated with EB can be explained within this 
DS model |^ 0, ^|. The importance of defects for 
the EB effect is also confirmed by recent experiments 
on Fe.r Zni .tF2 / Co bilayers and by investigations 
[T^ IT3 . Il5| where it was shown that it is possible to 
modify EB by means of irradiating an FeNi/FeMn sys- 
tem by He ions in presence of a magnetic field. On the 
other hand, models which assume the formation of a do- 
main wall parallel to the interface during field cooling 
cannot explain these findings. A domain wall parallel 
to the interface stores the energy for redirecting the FM 
magnetization like a spring that is wound up and non- 
magnetic defects will rather suppress this spring effect 
leading to a decrease of EB with increasing defect con- 
centration. Further support for the relevance of domains 
in EB systems is given by a direct spectroscopic observa- 
tion of AFM domains [TallTj. 

An important outstanding problem in EB systems is 
an explanation of the asymmetry of the reversal mode 
during hysteresis observed in some EB systems, either 
directly as an asymmetry of the shape of the hysteresis 
loop or - more detailed - by neutron scattering tech- 
niques |l8l Hgf . Even though in recent simulations an 
asymmetry has been observed a systematic theo- 
retical investigation of this effect is still lacking. In the 
following we will show that the occurrence of an asym- 
metry depends on the angle between magnetic field and 
easy axis of the AFM. A systematical variation of this an- 
gle reveals a rich variety of different reversal modes not 
seen so far in simulations or in theoretical investigations. 
These findings are not only important for a deeper un- 
derstanding of the EB mechanism but also essential for 
a correct interpretation of the complex reversal behavior 
observed experimentally. 

We consider the so-called DS model which was intro- 
duced recently, for a detailed discussion of its properties 
see In the following we have one FM monolayer 

exchange coupled to a diluted AFM film consisting of 
3 monolayers with dilution pmt — 0.5 for the interface 
layer and Pvoi = 0.6 for the remaining two layers (The 
geometry is sketched in ^3)- We deliberately chose these 
values since these yielded the largest EB in previous sim- 
ulations. 
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The FM is described by a classical Heisenberg model 
(spin variable S_^) with nearest neighbor exchange con- 
stant JpM- We introduce an easy axis in the FM (z-axis, 
anisotropy constant dz = 0.02Jfm)- The dipolar inter- 
action is approximated by an additional anisotropy term 
(anisotropy constant — — 0.2Jfm) which includes the 
shape anisotropy, leading to a magnetization which is 
preferentially in the y — z-plane. The AFM is modeled 
as a magnetically diluted Ising system with an easy axis 
parallel to that of the FM (Ising spin variable cr^ = ±1, 

= 0, 1 depending on whether site i carries a magnetic 
moment or not). Thus the Hamiltonian of our system is 
given by 

'H~-Jafm ^ eie-jUiU-j- ^ B^eiUi 

{i,i>6AFM isAFM 

- ^FM ^ ■ S_j - ^ {d^Siz + dxSf^ +R - 

(ij)eFM ieFM 

— Jmr ^ ^if^iSjz- 

(leAFMj'eFM) 

The first line describes the diluted AFM, the next line 
contains the energy contribution of the FM and the last 
line includes the exchange coupling across the interface 
between FM and AFM, where it is assumed that the 
Ising spins in the topmost layer of the AFM interact 
with the z component of the Heisenberg spins of the FM. 
For the nearest-neighbor exchange constant Jafm of the 
AFM which mainly determines its Neel temperature we 
set Jafm — — Jfm / 2 and we assume Jint — — Jafm ■ We 
use Monte Carlo methods with a heat-bath algorithm and 
single-spin flip methods for the simulation of the model 
above. The trial step of the spin update is a small vari- 
ation around the initial spin for the Heisenberg model 
and - as usual - a spin flip for the Ising model j2l|. We 
perform typically 136000 Monte Carlo steps per spin for 
a complete hysteresis loop, a number which turns out 
to be sufficient for the observed reversal modes to be in 
quasi-equilibrium. The lateral extensions of our system 
are Ly — Lz = 128 and periodical boundary conditions 
are used within the film plane. The external field S is 
within the film plane where 9 denotes its angle to the 
easy axis of the FM and AFM. 

With the FM initially magnetized along the z axis we 
first cool the system in an external field B = JpM with 
= 0° from k^T = Jfm to O.IJfm which is from above 
to below the ordering temperature of the AFM. Then, for 
each angle 6 we reset the FM spins to be aligned along 
the direction of the external applied field B and let the 
system relax. The initial conditions for different angles 
are identical, i. e., we use the same defect realization and 
cooling field orientation for each series of angles and ro- 
tate the applied field B after the initial field cooling pro- 
cess. The simulation of each hysteresis loop starts with a 
field B = 0.4Jfm and is reduced to —0.4 Jfm (decreasing 
branch) in steps of AB — 0.004Jfm before being raised 
again to its initial value (increasing branch). Hysteresis 
loops are simulated for angles between 0° and 72° with 



an increment of 4° each time taking a configurational 
average over six different defect realizations. 

Typical hysteresis loops for three different angles are 
depicted in Fig. ^ Shown is the projection m|| of the 
magnetization along the applied field _B (upper graph) as 
well as the perpendicular component m± (center). The 
lower graph shows the magnetization of the AFM inter- 
face layer (along z direction). From the perpendicular 
component one can already see the strong influence of 
6 on the reversal behavior. For small angles we get an 
EB that compares to the ones of previous simulations. 
With increasing angles however it becomes smaller and 
even turns positive. An explanation for this effect will be 
given later. First we will focus on the reversal modes. 
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FIG. 1: Projection of the magnetization along (upper graph) 
and perpendicular (central graph) to the applied field _B. The 
lower graph shows the AFM interface magnetization (along z 
direction) . 

Depending on the angle 9 one observes either a sym- 
metric reversal by coherent rotation (8°), an asymmetric 
reversal by coherent rotation for the descending branch 
and a non-uniform mode with zero perpendicular mag- 
netization for the ascending branch (60°), or symmetric 
reversal by coherent rotation (72°) but - in contrast to 
the usual behavior - with the reversal along the same 
direction for the descending and ascending branch. 

These effects are even more clearly shown in Fig. 12 
where the same angular dependence of the magnetiza- 
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tion reversal is plotted as versus ruy, with rUz {rriy) 
denoting the magnetization parallel (perpendicular) to 
the Ising axis of AFM and FM. Each plot starts at point 
1 corresponding to a magnetization with the maximum 
applied external field _B. Here the spins of the FM are 
aligned with the external field. With decreasing field _B 
the magnetization starts to rotate clockwise and the par- 
allel magnetization will eventually vanish thus yielding 
the coercive field B~ (point 2). For an increasing field 
the corresponding field value is i?+ (point 4). When a 
negative maximum field is applied (point 3) the spins are 
aligned antiparallel with respect to their initial orienta- 
tion at point 1. 
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FIG. 2: Magnetization my versus niz for different angles 8; 1 
and 3 indicate points with maximum external field \B_\ and 2 
and 4 points where m|| vanishes. The horizontal dashed line 
corresponds to the easy axis of FM and AFM and the other 
dashed line to the field axis. 

For 6 = 8° (Fig. [2^) we observe mainly coherent rota- 
tion for magnetization reversal on both sides of the hys- 
teresis loop. When the external field is either decreased 
or increased the FM spins first rotate towards their clos- 
est easy axis before reversal eventually occurs. Therefore 
one obtains different signs for the perpendicular magne- 
tization for decreasing and increasing fields, respectively. 
In Fig. I^Jd (0 = 60°) one observes that for a decreasing 
external field the system behaves as before displaying co- 
herent rotation. When the external field is increased, the 
system shows an entirely different behavior, though. The 



contribution of the magnetization perpendicular to the 
field axis is negligible. Obviously, here we have a highly 
non-uniform reversal mode with vanishing total magneti- 
zation. Increasing 9 even further to 72° (Fig.[2t) reveals 
again reversal mainly by coherent rotation on either side 
of the hysteresis loop with a smaller net contribution mj_ 
on the increasing branch. Most notably however is the 
fact that now the rotation occurs via the same side for 
both branches of the hysteresis loop, in contrast to the 
reversal for small angles. 



field I 



easy axis 




FIG. 3: Sketch of the effective field contributions from ex- 
ternal field, Hz, exchange field of AFM Hx, and uniaxial 
anisotropy Ha acting on the FM close to the coercive fields 
and B~ before the reversal sets in. 

For an understanding of these effects we note, that 
there are three different contributions to the mean ef- 
fective field acting on the FM during its rever- 
sal, namely the exchange field of the AFM, H ^ = 
Jint^t-afmI aligned with the easy axis of the AFM, the 
external magnetic field, -ff^ — and the anisotropy 
field 4 = 2D zm-pMZ. aligned with the easy axis of the 
FM, here, also the z axis. 

For the particular case shown in Fig. \Bp (60°) these 
effective field contributions are sketched in Fig.|3]for two 
characteristic values of the magnetization, a and b in Fig. 
12). As is shown in Fig.^), for this value of 9 the mag- 
netization curve of the AFM which is shifted upwards 
due to the fact that the AFM is in a DS with a surplus 
magnetization after field cooling 10] leads always to a 
positive effective exchange field H ^ . Additionally an ef- 
fective field coming from the anisotropy of the FM acts 
on the FM which tries to rotate the magnetization to its 
closest easy-axis direction. It depends on the sign of the 
FM magnetization and, consequently, points into differ- 
ent direction on both sides of the hysteresis. The fields 
labeled with — in Fig. |31correspond to point a in Fig. I^b, 
the fields with + to point b. In the first case the effective 
field has a large angle with the magnetization leading to 
a strong torque which favors coherent rotation. In the 
second case the effective field is more aligned with the 
magnetization which favors non-uniform reversal modes. 
Note that this is in particular true for 9 = 0° where the 
effective field is aligned with the applied field on both 
sides of the loop. Here, in agreement with the above 
discussion we observed in the simulations non-uniform 
reversal modes on both sides. However, we do not show 
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these results here, since we beheve that the perfect ahgn- 
ment of field and easy axis is a very special case which 
hardly will occur in an experimental situation. 

The strength of the anisotropy field depends on the 
projection of mpM onto the easy axis and decreases with 
increasing angle 9. When H\ is smaller than the 
reversal on the ascending loop branch will be over the 
same side as that of the descending branch (as in Fig. 
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FIG. 4: Coercive field values and B~ versus angle between 
field and easy axis of FM and AFM. 

A finite angle 9 between the applied magnetic field and 
the easy axis of the AFM can lead to another surprising 
effect: the possibility of a change of the sign of the bias 
field when changing this angle. The reason is a different 
dependence of the coercive fields on 9 due to the different 
reversal mechanism discussed above. The coercive fields 
as function of 9 are shown in Fig.^J For the reversal 
mechanism is always a coherent rotation and its value 
decreases with 9 monotonically. On the other hand, for 
_B+ the behavior is more difficult. Here we found a non- 
uniform reversal mode for intermediate angles for which 
increases with 9 which may even lead to change of 
sign of the bias field for certain angles. Note, however, 
that this requires bias fields which are rather small as 



compared to the coercive fields, so that this effect might 
not always be observed. 

To summarize, varying the angle between the applied 
magnetic field and the easy axis of FM and AFM we 
found either identical or different reversal modes on the 
ascending and descending branch of the hysteresis loop. 
For 9 = 60° the reversal is maximal asymmetric, with a 
coherent rotation for the descending branch and a non- 
uniform reversal with zero perpendicular magnetization 
for the ascending branch. We note, that different angular 
dependencies are obtained for different system parame- 
ters, e. g. varying the exchange constants or the AFM 
layer thicknesses. For a comparison with experimental 
situations one should note that in realistic systems the 
situation might be even more complicated due to i) a 
twinned or granular structure of the AFM leading to a 
less well defined anisotropy axis ii) a possible additional 
angle between easy axis of FM and AFM which in our 
investigation are identical. However, our investigations 
show that in general one has to deal with different re- 
versal modes for the ascending and descending branch, 
respectively, resulting in interesting and even surpris- 
ing new effects. Furthermore, our model reveals that a 
twinned crystal structure of the AFM is not necessary for 
the occurrence of an asymmetry in contrast to what was 
conjectured earlier p3 |. In corresponding experimental 
investigations a systematic variation of the angle between 
the magnetic field and the easy axis of FM and AFM is 
highly desirable. We expect a deeper insight into the 
physics of EB from a comparison of those results and our 
simulations. 
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